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SUMMARY 


The seismic analysis of the FFTF-PIOTA (Fast Flux Test Facility-Postirrad- 
ation Open Test Assembly), subjected to a horizontal DBE (Design Base Earth- 
uake) is presented. The PIOTA is the first in a set of open test assemblies 
•o be designed for the FFTF. Employing the direct method of transient analysis, 
he governing differential equations describing the motion of the system are set 
ip directly and are implicitly integrated numerically in time. A simple lumped- 
lass beam model of the FFTF which includes small clearances between core com- 
ponents is used as a "driver" for a fine mesh model of the PIOTA. The nonlinear 
-orces due to the impact of the core components and their effect on the PIOTA 
ire computed. 


INTRODUCTION 


The reactor core of the FFTF**is designed to accommodate bowing of fuel 
issemblies which is caused by the core's neutronic and thermal environment. 

'he individual fuel assemblies have a floating collar design at the above core 
load pad where adjacent assemblies contact to take into account the deforma- 
:ions of the core components. During certain periods of the reactor cycle 
.mail clearances may exist between core components; hence, it is important to 
ietermine the effect of design impact loads between these structural components 
to preclude unacceptable stress levels and failures. 

The topic of impact is the subject of numerous studies. Special purpose 
computer program solutions have been developed for analysis of impact of reactor 
internal structures by Bohm and Nahavandi , reference 1, using explicit inte- 
gration procedures. The analysis reported in this paper uses implicit inte- 
gration procedures through the transfer function capabilities of the general 
Durpose computer program NASTRAN to solve for the nonlinear loads due to impact 
between structural members. 

*The Hanford Engineering Development Laboratory is a United States Energy 
Research and Development Administration Laboratory. HEDL is operated bv 
the Westinghouse Hanford Company. 

**See Nomenclature Table. 
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NOMENCLATURE* 


FFTF-PIOTA 

DBE 

NASTRAN 

OTA 

DOF 

BANDIT 

PWR 

IVHM 

HEDL 

CYBER 


Fast Flux Test Facility-Postirradiation 
Open Test Assembly 

Design Base Earthquake 

NASA STRUCTURAL ANALYSIS 

Open Test Assembly 

Degrees of Freedom 

Computer Program — to determine minimum 
bandwidth 

Pressurized Water Reactor 

In-Vessel Fuel Handling Machine 

Hanford Engineering Development 
Laboratory 

HEDL Computer-Control Data Corporation 
Model 74-18 


TRD 

DMAP 

PNDL 

MSC 

INFONET 

COSMIC 

I/O 

CPU 

SDRC 

CYBERNET 


Transient Response Displacement 
Direct Matrix Abstraction Program 
A DMAP Module 

McNeal Schwendler Corporation 

Computer Sciences Corporation Computer 
Network System 

Agency for United States Government Release 
of Computer Programs 

Input/Output Units 

Central Processing Units 

Structural Dynamics Research Corporation 

Control Data Corporation Computer Network System 


*In order of jse. 
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SYMBOLS 


Jaa 

damping matrix 

'aa 

stiffness matrix 

iaa 

mass matrix 

ila 

1 

i 

a set of coordinates expressed as a column 
matrix of two terms for each grid point; 
y, e z (see figure 3) 

| 3 

pivotal frequency 

fit) 

seismic forcing function 

H n) 

nonlinear impact load 

9 

percent critical damping 


MODELING PROCEDURES 


Figures 1 and 2 portray the salient points of the FFTF reactor. The 
reactor is idealized with lumped-mass, spring, gap and beam elements for trans- 
lation and rotational motion. The DBE , a translational input is applied dir- 
ectly to the reactor head. The core components, fuel assemblies, etc., are 
pinned at the core barrel and lateral support is provided at the pad locations 
^ th6 j 9ap ?' ements * The core components, if permitted to rotate, may be con- 
sidered as inverted pendulums. The core components (reflector, fuel assembly 
3TA), figure 3, are expanded horizontally until a return is made to the core 
oarrel , for an enclosed representation of the FFTF core. 


Figures 3 and 4 show the grid point and finite element notation for the 
pourse FFTF reactor mpdel and the fine mesh PIOTA model. A fine mesh was used 
to model the PIOTA, figure 4, because of the boundary conditions between the 
PIOTA's components, the close coupling of its natural frequencies, and the 
predicted critical stress area at the PIOTA's outlet ports. The coarse reactor 
model is numbered in the 1000 series of grid points and 2000 for the element 
numbers. The model parameters such as size, mass, stiffness, and damping values 
are described for the respective models in references 2 and 3. The coarse 
reactor model is like the FFTF systems model used for both SCRAM and nonlinear 
seismic analysis, reference 2, with the exception that certain reactor 
components are not included in the subject model; i.e.. Instrument Tree, 

IVHM, etc. The SCRAM and nonlinear seismic analysis, reference 2, was per- 
formed using proprietary (Westinghouse Advanced Reactor Division) special 
purpose computer programs. 


The set of equations used to determine the seismic response is 
[Maa] {Ua} + [Baa] {Ua} + [Kaa] {Ua> = {f( t) } + {f (n) } 
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The physical set of equations for the models shown in figures 3 and 4 has 
927 DOF, of which 537 belong to the PIOTA. The DOF after reduction for 
coordinates y and e z , boundary conditions, constraints, etc., J, s a 2 f; ' 

problem size (297 DOF, time step increment of 0.0005 second and a 5 second 
earthquake) involves a sizeable computation investment, approximately 3 x 10 
iterations. Hence, further reduction of the problem was desirable. Of the ^ 
three reduction methods considered, substructuring, modal synthesis and Guy an, , 
the Guyan method was found best for minimum loss of accuracy in conjunction 
with nonlinear dynamic analysis. Its basis is that fewer D °F ar ® ne ®ded to 
describe the inertia of a structure than are needed to describe its elasticity 
with comparable accuracy. In using the Guyan method NASTRAN recovers directly, 
displacements, accelerations, etc., for the coordinates used in the reduction | 

process. The Guyan reduction method was used to reduce the 297 DOF to an 

analysis set of 104 DOF. The Guyan reduction process increased the bandwidth . 

from 3 to 45; however, the run time was reduced by 65% based on CPU, see 

table I. The preprocessor BANDIT was used to reduce the bandwidth for the | 
PIOTA model, figure 4. BANDIT was not applied to the FFTF model, figure 3, 
because of the gap elements and the straightforwardness of the model. 


NONLINEAR TRANSIENT ANALYSIS 


In this transient solution, the coupled equations of motion are inte- 
grated directly without any uncoupling by modal methods. Rigid Format No. 9. 

The two basic methods used for direct numerical integration are explicit 
integration in time and implicit integration in time. Difference formulas 
that relate the accelerations, velocities and displacements are used in both 
the explicit and implicit integration methods. NASTRAN (reference 4) uses a 
form of the Newmark Beta Method, implicit integration in time, that yields an 
unconditionally stable solution for a wide range of transient dynamic problems. 
The stability limit is a function of the period of the highest vibration mode 
of the system. Though the implicit integration method is not as fast per 
time step as the explicit method, the unconditional stability permits the use 
of large time steps. A time step of 0.0005 second was selected for the runs 
described herein and no numerical stability problem was encountered. The 
introduction of nonlinearities in the implicit method of integration may cause 
numerical stability problems in addition to those mentioned above. These prob- 
lems are due to the inconsistent definition of displacement and velocity 
between linear and nonlinear forces and may result in the presence of a para- 
sitic mass on the coordinates to which nonlinear forces are applied. The 
remedy is to add sufficient mass to the coordinate directly or to reduce the 
time step of integration so that the parasitic mass effect is negligible in 
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[omparison to the nodal mass. A time step of 0.005 second was used 

to obtain a linear solution for the model shown in figure 3. 

r his is the same time interval as the input acceleration -time history. With 

:he addition of gap elements to the model, figure 3, the time step was reduced 

po 0.0005 second through the above parasitic mass consideration. 

, In NASTRAN, the nonlinear effects (gap elements) are treated through the 
jse of an additional applied load vector. The gap element and the respective 
ipplied load vector are user-created by means of transfer functions. Gap 
ilements similar to the ones shown in figure 3 are described in reference 5 
nth transfer function and nonlinear load card images. 


DAMPING AND THE GAP ELEMENT 


Two forms of damping* are used in this seismic analysis, structural and 
impact. In both forms, assumptions concerning the effects of damping on the 
lodal coordinates are based on a viscous model in which the energy dissipated 
per cycle is proportional to frequency. A uniform structural damping of 2% 
of the critical for the DBE was input to the model (figure 2) in terms of 
Stiffness, as follows: 


[Baa] = 2^ 
Wo 


[Kaa] 


This method, reference 4, of inputting equivalent viscous damping is an approx- 
imation, since the viscous damping forces are larger at higher frequencies and 
smaller at lower frequencies. The structural damping. Baa, has small effect, 
or is of no effect, on the response of the model. This is due to the over- 
ruling effect and nature of impact damping (10% vs 2% of the critical). 


The impact damping is related to the nonlinearities of the gap element. 

The impact damping and stiffness, c and k, figure 5, are based on a coefficient 
restitution method (rebound) and the Hertz theory of impact of two solids, 
i.e., an elastic statical consideration. The relationship between coefficient 
of restitution and critical damping is derived in reference 2 and is shown in 
.this report as figure 6. An impact damping value of 10% of the critical was 
used in this analysis. This represents high viscous damping. Impact damping 
values for FFTF reflector and fuel assemblies have not been measured exper- 
imentally. However, there are data available from PWR fuel assemblies tests, 
reference 2, and from a FFTF IVHM test, reference 5, to corroborate the high 
impact damping value. Figure 6 shows a comparison of these impact damping 
values, 19% for PWR fuel assemblies and 15% for the vertical IVHM test. The 


*Structural - hysteretic damping 
Impact - Maxwell representation of viscous damping 
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IVHM 15% value was determined from the rate of decay of successive re- i 

bounds. Therefore, the 10% of the critical damping is considered conser- 
vative for these analyses. 

The vertical drop IVHM test, reference 5, affords a means to assess the 
validity of the NASTRAN gap element. In comparing the experimental results 
with those computed from the simple scalar lumped-mass NASTRAN model, there 
was close agreement. The free fall times and displacements for this test and 
model were readily calculated by hand methods. 

i 

The general form of the gap element used in this analysis is shown in 
figure 5(a). Gaps in positive and negative directions (for example, the 
boundary conditions of a rod moving within a tube), can be represented by 
this element. A gap element with closing capability only (i.e., that portion I 
of the force deflection curve shown in the third quadrant, figure 5(b), and the 
analogous mechanical model shown in figure 5(c)), was used to represent the 
clearances between the core components, at twelve places (figure 3). The ex- j 
panded representation of the core with the gap closure described above, yields j 
a closed core configuration for impact analysis of individual core components. 
Figure 5(c) shows the mechanical model where the impact spring and gap are j 

accompanied by energy dissipation, and viscous damping. Functionally, the 
viscous damper, Cx should act only when the gap is closed, x : >u x . At the : 
time these seismic analyses were performed, a switchable, viscous damper had 
not been verified for use in NASTRAN. A pseudo Maxwell model, figure 5(c), was 
used in this analysis; the Maxwell model does not necessitate a switchable 
damper, since the damper Ci is in series with the impact spring k x . The 
linear spring k 3 , which was small compared with k ls improved the stability of 
the numerical solution, i.e., a larger time step was used. 


SEISMIC TIME HISTORY 


The input to the FFTF seismic model. Grid Point 1001, figure 3, consisted 
of an acceleration time history with data points at 0.005 second intervals. 
This seismic transient is the singular output at the concrete ledge from a 
two-dimensional finite element soil interaction model of the containment 
building (reference 6). The dynamic coupling between the containment building 
and the reactor vessel is assumed to be negligible. 

A DBE with 2% of the critical damping, and a duration of 20 seconds was 
specified (reference 6). Previous investigations (reference 2) revealed that 
only the accelerations from 1.5 seconds to 6.5 seconds were significant; this 
then is the seismic transient that was used for the PI0TA analysis (figure 7). 
The DBE earthquake was divided into ten increments, as shown in figure 7, for 
computer restart advantage and clarity of output. The printout time interval 
was ten times the integration time step, or 0.005 seconds. The analysis was 
started at 1.5 seconds; thus, there is a 1.5 second time shift in the output 
response, see figure 7. 
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SOFTWARE AND COMPUTER FACILITIES COMPARISON 


This section compares the off-load* results with those obtained on the 
IEDL-CYBER for a 0.5 second DBE input to the FFTF-PIOTA model. The CYBER 
fas not available to perform the 5 second DBE due to other commitments. 

In performing the CYBER and off-load benchmark runs we found that the 
)0SMIC version of NASTRAN, Level 15.5, computed incorrect nonlinear loads. 

[he nonlinear load computations were evaluated by the IVHM model described 
n reference 5; the results for the IVHM model can readily be checked by 
tand methods. The IVHM evaluation revealed that the PNLD module in NASTRAN 
Level 15.5 was not functioning correctly. DMAP instruction number 139, 
IASTRAN Rigid Format No. 9 gives the relationship of the PNLD module to the 
imputations of TRD. The source language for NASTRAN is essentially written 
n FORTRAN IV; the DMAP instructions are a compilation of this source 
anguage. The IVHM analysis was successfully performed in 1972 on the 
INIVAC 1108 using Level 15.1, hence a CDC 6600 version of NASTRAN Level 15.1 
vas installed on the CYBER and used for the CYBER computations described in 
this report. 

Benchmark runs were performed at the four computing facilities shown on 
table I. The CYBER run was used as a reference for evaluating the off-load 
results . The reference run column II, table I, is identical to the FFTF- 
5 I0TA runs described herein with the exception of the OMIT 1 cards. The 
3MIT bulk data cards are used to achieve the Guyan reduction described 
oreviously in this report. It is difficult to compare results from dis- 
similar computing facilities because of timing algorithms, input-output 
variances, system dependent software, etc. These problems were alleviated 
in the FFTF-PIOTA evaluation by comparing the run time required for the TRD 
module, see table I; this module uses the majority of the central processing 
time needed for direct transient analyses. 

Table I, columns III through VI, show that the MSC version of NASTRAN 
gives significant reductions in run time when compared to the HEDL or INFONET 
COSMIC version of NASTRAN. The benchmark run that the MSC performed in 
Los Angeles, column IV, had two alterations, the transfer functions used for 
the gap elements were replaced by multipoint constraint equations to utilize 
the more efficient symmetric matrix decomposition rather than unsymmetric 
decomposition, and certain data that was transferred to peripheral storage in 
Level 15.1 was held in main core in MSC-20, thus eliminating peripheral 
[processing-calls. The above alterations account for the reduction in I/O 
time of 1237 to 26 seconds and CPU reduction of 864 to 569 seconds. The MSC- 
Los Angeles results shown in column IV were obtained subsequent to the 
analyses described in this report. 


*0ff-load — Computer facilities other than available at HEDL-CYBER 
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The deck used by SDRC in their benchmark run was identical to the CYBER 
run, columns II and III, table I. The improvement in run time, I/O 3300 
versus 1237 seconds and CPU 1134 versus 864 seconds, is due primarily to the 
version of NASTRAN utilized, COSMIC Level 15.1 versus MSC-20; the operating 
systems, CYBER and CYBERNET, are similar, both are basically CDC 6600 machines. 
The solution for the SDRC benchmark run was identical to the CYBER results, 
columns II and III. Therefore, as a result of these benchmark evaluations 
and the need for off-load capability, the SDRC facilities at Cincinnati, Ohio, 
NASTRAN MSC-20 and the CYBERNET operating system were used to perform the 
5 second earthquake run described in this report, see column VI, table I. 


RESULTS 


NASTRAN, as a large, general-purpose, structural -analysis program, has 
features, through checkpoint and restart, to recover and compute element 
stresses, forces and moments for each of the elements shown in figures 3 
and 4. Output for this analysis was restricted to bending moment and shears 
for the boundary-adjacent elements numbers 1, 24 and 52, and for element 
number 46, a reduced section on the output port of the holddown tube, 
figure 4. Tabulated peak values for each increment are shown in table II. 
Figures 8 through 11 show typical load versus time plots. The amplitude of 
these PIOTA loads are largely dependent on the nonlinear loads caused by 
impact between core components, see figure 9. The effect of the impact 
load on the wave form of the element loads can be seen by noting that core 
impacting at Grid Point 1022 does not occur until 0.18 seconds, figure 9; 
the corresponding loads in Element 52, figure 11, are very low until this 
time. The wave form of the PIOTA loads at the head, figure 10, is not 
directly influenced by the impacting of core components, the PIOTA acts as a 
mechanical filter. 


CONCLUDING REMARKS 


A seismic analysis of a reactor internal component with impact between 
core components has been performed. From this effort, it can be concluded 
that NASTRAN has good nonlinear transient analysis capabilities. At the 
present time, solutions are relatively expensive from a computational stand- 
point when compared with solutions from special purpose computer programs. 
However, with the efficiencies projected for Level 16 and the use of certain 
operating system and modeling procedures noted in this report, NASTRAN can 
effectively be used for nonlinear seismic analysis. 

In using NASTRAN's nonlinear features, it is suggested that one begin 
with small sample problems with known solutions. The setup of nonlinear 
elements using NASTRAN's transfer function capabilities is a user-oriented 
function which presents opportunities for data errors; data checking by 
NASTRAN is minimal since the nonlinear features are mathematical abstractions. 

The design of the PIOTA was found to be structurally adequate for the DBE. 
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FFTF-PIOTA FINITE ELEMENT MODEL— SEISMIC BENCHMARK- COMPARISON — NASTRAN RUN TIME-MODULE TRD 
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SUMMARY-PIOTA 
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MAXIMUM VALUES FOR EACH INCREMENT, BASED ON M J 
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Figure 2.- Schematic of FFTF lumped-mass beam model. 
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Figure 5.- Gap element. 
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8.- Increment 1 - acceleration input (see figure 7). 
1 in/sec 2 = 0.0254 m/sec 2 . 
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Figure 10.- Increment 1 - moments, element 1 - PIOTA 
(see figure 4). 1 in-lb = 0.11298 N-m. 




